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Two-particle like-sign and unlike-sign correlations including Coulomb final state 
interactions are calculated for Gaussian wave packets emitted from a Gaussian 
source. We show that the width of the wave packets can be fully absorbed into 
the spatial and momentum space widths of an effective emission function for plane 
wave states, and that Coulomb final state interaction effects are sensitive only 
to the latter, but not to the wave packet width itself. Results from analytical 
and numerical calculations are compared with recently published work by other 
authors. 



1. Introduction. - To analyze the geometry and dynamics of the collision 
region, two-particle correlations C(q, K) of like-sign and unlike-sign hadrons 
have been studied extensively in relativistic heavy ion collisions at AGS 111 
and CERN SPS HH energies. They show q-dependent structures at relative 
pair momenta \q\ < 100 MeV. These originate mainly □ from (i) final state 
interactions (which for pions at \q\ < 100 MeV are dominated completely by 
the Coulomb force) and (ii) the quantum statistics of identical particles. 

Practical attempts to reconstruct space-time information from two-pion 
correlation data in momentum space so far exploit mostly the quantum sta- 
tistical Hanbury Brown - Twiss (HBT) effect between identical bosons!!!!!]!. 
This requires a prior subtraction of final state Coulomb interaction effects from 
the measured correlation functions, with proper account for the finite size of 
the emission region @0§ illl3O. Up to now this is done directly in the experi- 
mental analysis, either by taking experimental unlike-sign pion correlations to 
correct the like-sign onesll, or by a model calculation for the Coulomb effect 
expected for a finite size emission region 0. 

This approach was recently questioned by Merlitz and PelteEHIll'lll. From 
a numerical analysis based on Gaussian wave packets emitted from a Gaus- 
sian source, they concluded thatllH "the expected Coulomb distortion in the 
momentum correlation . . . becomes unobservable" and that therefore "exper- 
imental data, which are published after Coulomb correction, are wrong for 
small momentum differences" . If correct, this conclusion would invalidate a 
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substantial part of the existing work on the analysis of two-particle correlation 
data since it implies that either (i) any attempt to base a space-time inter- 
pretation of identical two-particle correlations on Coulomb corrected data is 
ill-founded or that (ii) any attempt to describe the particle emitting source in 
heavy ion collisions by a set of Gaussian wave packets is inconsistent. 

This dramatic perspective has led us to reconsider the calculation of two- 
particle correlations for Gaussian wave packets. Following Merlitz and Pelte we 
describe the particle emitting source in terms of a distribution of wave packet 
centers and a characteristic wave packet width a of the emitted particles. The 
emitted Gaussian wave packets are propagated into the detector under the 
influence of mutual Coulomb final state interactions. We derive analytical ex- 
pressions which show that the wave packet size a can always be absorbed by a 
redefinition of the model parameters characterizing the source size. For Gaus- 
sian source models, two-particle correlation measurements cannot differentiate 
between "source size" and "wave packet width" . In a limiting case we further 
prove analytically the equivalence of the Gaussian wave packet formalism and 
the usually adopted plane wave calculations irrespective of the size of the wave 
packet width. These analytical calculations show quite generally that the prob- 
lem pointed out by Merlitz and Pelte does not exist. Their results disagree 
with our numerical calculations as well as with analytical formulae which we 

derive without apnroximations from the same starting point as the calculation 
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presented in Ref.lif . 



2. Unlike-sign pion mrrelaiinn s. - Our starting point is a set of Gaussian 
one-particle wave packets Ii^'t3c3li3 

f (t tn) - - e -(3>-r i f/(2(T 2 )+ip i -se (i\ 

Ji{X,l ) - ^ cr 2)3/4 e ' W 

which are centered at initial time t = to at phase-space points (f*i,p,J. We ex- 
pand the time evolution of the corresponding two-particle state (xx, X2, to) 
= fi(xx,to) fj(x2,to) in terms of plane waves 4> Pl , P2 , 

J^y Aj{Pl,P2,t)<p PuP2 {x 1 ,X 2 ,t) , (2) 

cf >PllP2 (x 1 ,X2,t) = e- lEt ^2 K (X) ( f> q/ 2(r)^e-^ +E ^e 2lK - x e^ r ,(3) 

which we write in terms of center of mass coordinates X = \{x\ + x%), K = 
hiPi + P2), an d relative coordinates r = {x\ — X2), q = {pi — P2)- The 
probability T-y for detecting at time t — > 00 the two particles prepared in the 
state \Py with momenta pi and P2 is given byE 

VnipuPl) - P ij (q,K)= lim \A* 3 ( Pl ,P2,t)\\ (4) 
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t — >00 >OG 



*?n<«+^/2l*y>- (5) 



Here we separated the state (xi, X2, to) into relative and center of mass 
wave functions, 

*ifw = -4^ e - (x -^ )2/CT2+2ii ^- x , (6) 

V > f %a 2\3/4 K ' \'J 

where Xij , Kij , , are the corresponding center-of-mass and relative co- 
ordinates constructed from the wave packet centers. In (||) we have also intro- 
duced the M0ller scattering operator Q + = lim^oo e lff (*-*<>) e -tH (t-t ) f Qr ^ e 
final state interaction hamiltonian H = Ho+V(r), Hq = -^Ax-^A r . This 
M0ller operator maps the plane wave 4> q /2 onto the solution of the Lippmann- 
Schwinger equation for the corresponding stationary scattering problem. For 
two-particle Coulomb interactions this is the Coulomb scattering wave 

(fi+^ /2 )(r) - ^/i(r)=T(l-ir,)e^ei"- r F(ir,;l;iz-), (8) 

2 

z± = \{qr±q-r), r}=—, (9) 

where e 2 /47r = a = 1/137, r — \r\, q = \q\, F(irj;l;iz-) is the confluent 
hypergeometric function, and 77 is the Sommerfeld parameter. Eq. (|^) applies 
for pairs with opposite charges; for like-sign pairs one replaces n 1— » —77. With 
the help of Eq. (||), the calculation of the amplitudes (|) is reduced to a six- 
dimensional integral. 

2.1. Gaussian source model. - We consider a toy model of simultaneous 
particle emission at time t = to for which initially the wave packet centers are 
distributed with Gaussians of widths R and A in coordinate and momentum 
space, respectively. It can be specified by the following normalized distributions 
of relative distances and pair coordinates: 

^(r,q) = ¥ ^ E e-^-^, S pail (X,k) = J -^e-^. (10) 
With the choice R 2 — R 2 /2, A 2 = mT, and a 2 = 2a 2 , this model coincides 



1 9 

with the one considered by Merlitz and Peltelifl We calculate the unlike-sign 



two-particle correlator via the two-particle spectrum Vij(pi,p2) of Eq. (0), 
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averaged over the distributions (|T^) and normalized to the corresponding spec- 
trum for pairs of non- interacting particles: 

C+-(gK)- pnt ^ K) (11) 

I(q,K) = J d% j d 3 q ij d 3 X ij d 3 K ij S le i(f ij ,p ij )S paiT (X i ^ 

(12) 

For the non-interacting case Eq. ( |l2] ) must be evaluated with V^° nint which is 
obtained by replacing in the amplitude (^J) the relative Coulomb wave f2 + 4> q /2 
by the plane wave 4> q /2- Since the center of mass coordinate is not affected by 
the two-particle final state interaction, the integrations over <Sp a ir(-^y ■, Ki 3 ■) 
drop out in the ratio ([ll]), and the correlator C H (q, K) becomes independent 
of the pair momentum K. One finds 

C+ ~^) =G ^)(-^) e ^ J <Pre-f^ + i rq F{ir,-l;iz^) 

d 3 r' e-£fp-i r '- q [F(ir);l;iz'_)}* e'^~^) {r ~ r ' )2 , (13) 



where G{rj) = \T{1 + ir))e i m i^ — 2iri] / (e 2 * n — 1) is the Gamow factor. It is 
important that this correlator depends only on the parameter combinations 

R 2 = R 2 + ^ , A 2 = A 2 + ^ . (14) 

This shows that for Gaussian models of particle emission, the wave packet 
width er can be absorbed in a redefinition of the model parameters. There is 
no measurement which allows to determine a independent of R and A. Of 
course, the specific c-dependence in ( |l4| ) still constrains the values which R 
and A can take. In particular, R, A always satisfy the uncertainty relation 
RA > Tl/2, rendering the exponent of the last term in (|l3| ) always negative. 

In the absence of final state interactions, the dependence of the correlator 
on the parameter combinations (|l4| ) has been noted repea tedlyHMlM The 
momentum spect ra and correlations are entirely determined by the "effective" 
emission functionll3't2l 

Sas(r,p) = J d 3 fd 3 pS(f,p)S w . P Xr -f,p-p), (15) 

which is a folding integral between the distribution of wavepacket centers 
S(f,p) and the Wigner density S w . p . of a single particle wave packet. For 
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Gaussian source models, S e fi depends on R and A onlyH^!. The same holds 
true in the presence of final state interactions where the correlator can be 
written for arbitrary model distributions Si r, p ) as a quite involved expression 
depending on S cff only (see Eq. (60) in Ref.lHI). Our result ( |i~3| ) is an explicit 
representation of this general relation, obtained for the Gaussian source models 
(|l0|). It allows further analytical and numerical studies: 

2.2. Limiting cases. - In two interesting limits the correlator C" 1 (q) can 
be further simplified analytically. Using hin^^oo (A 2 /47r) 3 / 2 exp ~^-(r — r') 2 
= S^(r-r') we find 



A lim oo C+ ~ (9) = WWW J l ^ {r)l2 ' (16) 

This expression, first written down by KoOPmEa jfi the .uaially adopted start- 
ing point for plane wave ca lcul ations EJflBH3'B& PJc Ht^rE^I; it was shown by 
Baym and Braun-Munzinger 113 to be well approximated by a semi-classical 
approach. The limit A — > oo can be taken for arbitrary values of the wave 
packet width a and is equivalent to the limit A — > oo which describes an emis- 
sion function S without momentum dependence. In this limit, the correlator 
(q) for simultaneously emitted Gaussian wave packets c oincide s exactly 
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with the starting point of conventional plane wave calculations irrespec- 
tive of the size a of the wave packet. Rescaling a then simply amounts to a 
change of the effective spatial size R of the source. 



What happens if the source size becomes large? Changing in Eq. (13) 
the integration variables r — > VRr, r' — ► V~Rr', and replacing the Coulomb 
wave function by its leading contribution, fVSrJ — * exp ^|Vi?r • q + 

ir) In (qr - q-r))^ +o(^),wefind 

lim C + -(q) = 1. (17) 

This is expected: as the source becomes larger, the average spatial separation 
between particles increases and their Coulomb attraction decreases, leading to 
a flat correlator in the limit of infinite source size. 

2.3. Numerical results. - One may wonder whether a large but realistic 
effective source size R can come sufficiently cl ose to the limiting case R — > oo of 
( |l7j ) to support the claim of Merlitz and PelteHSEL^I that the Coulomb repulsion 
becomes effectively unobservable. To study this question we have calculated 
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the correlator (|T^) numerically, after doing the azimuthal integrations (q = \q\): 



C + -(q) = 47r 2 G(-r?) 



( A 



4ttR 



'■ dr e si 3 



dxei qrx F(ii T ,l;^qr(l-x)) 



o J -i 

oo /2 />1 

r' 2 dr'e~i^ / dy e-^ qr ' y [F(ir]; 1; |gr'(l - y))}* 

i! 7-1 



x J, 



(2B 2 rr'VT^Vl^ Ij 2 ) e -B 2 (r 2 -2rr'x 1/ +r' 2 ) 



(18) 



Here J is the modified Bessel function and B 2 — A 2 /4 + 1/(16.R 2 ). The 
numerical results for (q) are shown in Fig. |l| for the model parameters 
R = 3.5 fm, A = 84 MeV and different values of the wave packet width a. For 
ct = 2.5 fm, these values correspond to the model parameters chosen in Ref.O. 

1 .5 r 
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R = 3.5fm, A =84 MeV 
calculation via eq. (19) - 
A -> oo limit with radius R - 
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Figure 1: Two-particle correlator of unlike-sign pion pairs for the Gaussian model jlo[). 
The correlation depends only on the model parameter combinations R 2 = R 2 + a 2 /2 and 
A 2 = A 2 + 1/(2<t 2 ). The results for A = 84 MeV agrees well with Koonin's expression ( [l6| ) 
obtained in the limit A — > oo. 

One sees that, with increasing source size R 2 = R 2 + a 2 /2, the resulting 
Coulomb correlation indeed becomes smaller than the Gamow factor, but it 
remains clearly observable even for a very large wave packet size a = 14 fm. We 
also compare in Fig. 0the full correlator C+-{q) of Eqs. (||/||) to the limit 
(|l6|), calculated for the same value of R. At least for the model parameters 
studied here, both expressions agree almost exactly. There seems to be no 
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possibility to distinguish on the basis of final state correlations between the 
Gaussian wave packet formalism and the plane wave calculation based on the 
approximate expression ([16]) which was used in other studies lilot3. As long 
as the Gaussian wave packet width a is included consistently in the definition 
of the source size R = \/R 2 + cr 2 /2, both formalisms lead qualitatively and 
quantitatively to the same result. We conclude that the differing results derived 
in Ref.ELHI with the help of a numerical simulation of the time evolution of wave 
packets are incorrect. 

3. Like-sign pion correlations. - For pairs of pions of identical charge 
the Coulomb final state effects are superimposed on the quantum statistical 
effects resulting from the symmetrization of the two-particle wave function. 
Paralleling the calculation of section 2. with Bose-Einstein symmetrized Gaus- 
sian wavepackets "J^ and plane waves 4> Pl , P2 , one arrives at the symmetrized 
asymptotic two-particle amplitude 

Urn Af*( Pl ,p 2 ,t) = (0 K I *r r ) [< n + I *« ) + I *y >] > (19) 

from which the two-particle momentum-space probability V^ E (pi,p 2 ) is again 
calculated according to Eq. (0). Averaging Vf E according to ( |i"2| ) over the 
model distribution gives the numerator of the two-particle correlator which 
we call 7 BE,Illt (q, K). We normalize it by the method of "mixed pairs": an 
uncorrelated (mixed) pair is described by an unsymmetrized product state 
without Coulomb interaction and leads, after averaging over the model distri- 
bution, to / nonint (qr, K) (see Eqs. Taking both distinguishable states 
fi{xi,t ) fj(x 2 ,t ) and fi(x 2 ,t ) fj(xi,t ) into account we have 

rBE,intf 

K) = /nonint(g; K)+ %1 { ^ R) = + C-(g) ■ (20) 



For the Gaussian model (10), the center of mass coordinate is affected nei- 
ther by two-particle final state interactions nor by two-particle Bose-Einstein 
symmetrization. The correlator ( p0| ) hence does not depend on the pair mo- 
mentum K. It splits into two contributions. The "direct term" can be obtained 
from C H (q) in Eq. ( [l3|) by changing the sign of the Sommerfeld parameter, 
— ?y — > rj. The "exchange term" is given by 

C cx {q) = Re (g(ti) (J^fi) J d 3 r e ~IW+i r -i F{-i m X;izJ) 

x /d 3 r'e~i^ + * r '- 9 [F(- 2 77;l;i2;)]*e"(^"^) (T - r ' )2 \. (21) 
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This integral can be simplified to a 4-dimensional expression similar to (fL8|). 
The limits A — * oo and R — ► oo of the first term C dlr (q) are obtained from 
by replacing — t] — > 7/. The corresponding limits for the exchange term 
are given by 



lim C ex (<?) 

A^oo 

lim C cx (<j) 

R^oo 



1 



(47Ti? 2 ) 3 / 2 



d 3 r 



cos(r • q) \^(r)\ 2 



= 0. 



(22) 
(23) 



As in the case of unlike-sign correlations, the correlator C ++ (q) depends only 
on the parameter combinations R and A, but not explicitly on the wave packet 
width er. In the limit A — > oo, the Gaussian wave packet formalism again 
coincides with the Koonin formula^ (now with the additional symmetrization 
factor 1 + cos(q • r) under the integral) which is the starting point of most 
conventional plane wave calculations. 




100 120 
q in MeV 

Figure 2: Two-particle correlator of like-sign pion pairs for the Gaussian model (|lo|). The 
correlation depends only on the model parameter combinations R 2 = R? + a 2 /2 and A 2 = 
A 2 + l/(2cr 2 ). The full correlator is well approximated by the standard Koonin expression 
for a static momentum-independent source of radius squared R 2 — 1/(4A 2 ). 

In Fig. H we compare numerically the full correlator C ++ (q) from Eq. (|o|) , 
for A = 84 MeV, with the limit A — ► oo, both calculated for the same value 
of R 2 = R 2 + er 2 /2. (Due to spherical symmetry of the source the correlator 
depends only on q = |q|.) For small values of a one observes a small, but 
significant difference. The reason is that even in the absence of Coulomb final 
state interactions, the HBT radius parameter (which gives the g-width of the 
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correlator) is not exactly given by the source size R 2 but rather by HI 

^HBT = i? 2 -^- (24) 

For large values of R or a, the term R 2 dominates this expression, and the 
difference between -Rhbt ana ^ disappears (see Fig. ||). In fact, when com- 
puting the limit A — > oo of C ++ (q) using -Rhbt instead of R 2 , the agreement 
with the full correlator C ++ (q) becomes almost exact even for small values of 
a (see inset in Fig. ||). We note that even for the smallest value studied here, 
a = 2.5 fm, the term 1/(4A 2 ) contributes only w 5% to -Rhbt- This small 
difference is clearly visible in the exchange term (|2l| ) of C ++ (q), whereas its 
influence on the direct term C dlr (and hence on the correlator C H ) is found 
numerically to be an order of magnitude smaller. This illustrates that the 
two-particle correlator C ++ of identical pions is more sensitive than C H to a 
small change in the Gaussian width of the phase space density. 

The modification ( pi| ) of the radius parameter to be used as input in the 
plane wave calculation can also be obtained- from the Koonin expression H 
or, most explicitly, from Eq. (65) of Ref. ELII. The consistency of Koonin's 
expression with the full correlator in the present model calculation is anon- 
trivial check of the so-called smoothness approximation used in Ref. liil to 
derive Koonin's expression from a general treatment of two-body final state 
interactions. 

To sum up: as long as the Gaussian wave packet width a is included 
consisterit.lv in the defini tion of the source size, both the plane wav e calcula- 

lS M4RPlJ4 |i0lpl| )22|,P3| QTlr | f>, riomooioTi were, r^W f™-™ai;o™ll6||l7|,|l8lfli 



tions liHflBLJa and the Gaussian wave packet formalism lil 

lead to qualitatively and quantitatively equivalent results. While the present 
study proved this only for Gaussian source models, we expect it to be true 
quite generally since we know that Ell the relation ( |l5|) between the effective 
emission function and the Wigner density of single particle wave packets holds 
for arbitrary model distributions and that 113 two-particle momentum corre- 
lations are mostly sensitive to the Gaussian characteristics of the source in 
space-time. 
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